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Abstract 

A planet orbiting around a pulsar would be immersed in an ultra-relativistic 
under-dense plasma flow. It would behave as a unipolar inductor, with a significant 
potential drop along the planet. As for Io in Jupiter's magnetosphere, there would 
be two stationary Alfven waves, the Alfven wings (AW), attached to the planet. 
The AW would be supported by strong electric currents, in some circumstances 
comparable to those of a pulsar. It would be a cause of powerful radio waves 
emitted all along the AW, and highly collimated through relativistic aberration. 
There would be a chance to detect these radio-emissions from Earth. The emission 
would be pulses as for ordinary pulsars; their occurrence would depend on the 
planet-star-observer angle. These results are still preliminary, further work needs 
to be done. 



1 Exoplanets around pulsars 

Wolszczan and Frail [1992] and Thorsett et al. [1993] discovered two planetary systems 
around a pulsar, through pulsar timing method (Doppler shift of the pulsar's period P). 
The data about the two stars, and their planetary systems are given in Tables [T] and [21 
Possibly, other planets around other pulsars remain to be discovered. In this paper, we 
investigate the ability of planets around pulsars to emit radio waves. 

Indeed, the flux of energy available in a pulsar's environment is huge, and even comparable 
to the gravitational energy that bounds the planet to its star. The pulsars are known to 
spin down, and the associated dissipated power is E rota tion = —MjQ^Q^ = Att 2 MiP^/P^, 
where is the neutron star's rotation velocity, P = 2^/0* the corresponding period, 
Q t . < is the rotation velocity time derivative, Mj ~ (2/5)M*i? 2 is the star's moment 
of inertia, M* is the neutron star's mass, and R* is its radius. In the case of PSR 
1257+12 that hosts three planets, E rotation = 2. x 10 27 W. The gravitational energy of 
planet "a" is Eq = GM*M p /2a = 4. x 10 32 J. Let us assume that this planet cap- 
ture the whole energy emitted from the pulsar through it section ttR 2 p , then the ratio 
(R p /a) 2 (E G /E rotation ) = 1,5 x 10 6 yr indicates that most of the gravitational energy of 
the planet would be dissipated in a million year. Of course, we do not expect that the 
energy capture process would be 100 % efficient, but this comparison proves that even 
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with a low efficiency, there may be enough power to trigger plasma acceleration around 
the planet, and powerful radio- waves emissions. 

There exists two main families of pulsars. The "standard pulsars" have a period P* ~ 1 
s, and P ~ 10~ 14 - 10~ 12 . The second family, to which belong PSR 1257+12 and PSR 
1620-26 are characterized by a period P* ~ 10 ms and P ~ 10~ 20 — 10~ 17 . This is the 
family of the "millisecond pulsars". In the following pages, we will refer to these two 
families of pulsars for numerical applications. 

A planet is immersed in the pulsar's magnetosphere, and the interaction between the 
pulsar's dissipated energy and the planet is supposed to be of electrodynamic nature. To 
understand it, we need to know in which region of the magnetosphere is the planet. The 
magnetosphere of a pulsar comprises two regions bordered by a cylinder, called the "light 
cylinder" that is parallel to the star's rotation axis, with a radius r^c = c/Q*, where c 
is the light velocity. It corresponds to an absolute limit : at the distance tlc, a plasma 
co- rotating with the star would have the velocity of light. Up to a few light cylinder radii, 
the structure of a pulsar magnetosphere is complex. Beyond, it is a wind, blown away 
from the inner magnetosphere, at relativistic velocities. 

Is a planet in the inner magnetosphere or in the wind ? Table [3] shows the radii and the 
corresponding keplerian orbital periods for the light cylinder of 1 s and a 10 ms pulsars. 
We can see that a planet orbiting inside the light cylinder would have an orbital period 
not exceeding 3 minutes, (being extremely close to the star) while the period of the known 
planets exceed several days. Therefore, it is reasonable to consider that planets orbit into 
the pulsar wind, far from the light cylinder. For instance, the distance of planet PSR 
1257+21 "a" is a = 9.4 x lO 4 ?^. If this planet was orbiting around a P = 1 s pulsar, the 
ratio a/r LC would still be large, of the order of 600. 



2 Alfven wings 

In [Mottez and Heyvaerts, 2011a], now notated [MH], the electromagnetic interaction of a 
relativistic stellar wind with a planet or a smaller body in orbit around the star is inves- 
tigated, through analytical calculations in MHD. Special relativity effects are included. 

For simplicity, the theory exposed in [MH] and summarised in this section, is developed 
for a pulsar with a magnetic dipole axis aligned with the rotation axis. The wind flow 
is supposed to be radial (r direction in spherical coordinates), and the magnetic field in 
the companion's environment is supposed to be azimuthal (cf> direction). In most wind's 
models for aligned pulsars (see [Kirk et al. 2009] for a review), the wind is characterized 
by two invariants along its flow: the neutron star's magnetic flux and the mass flux /, 
defined as 

/ = 7oP» 2 , (1) 
* = r 2 Bl (2) 

where 7 is the wind's Lorentz factor associated to the unperturbed wind's velocity v , and 
Bq is the radial magnetic field. The Lorentz factor is expected to be high, but estimates, 
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based on various observations and models, vary over a fairly large range from 10 1 to 10 7 
[Kirk et al. 2009] . The condition E + v x B = implies that the azimuthal magnetic field 
is given by, 

B* = B ?" - ^ ~ -go^V. ( 3 ) 

The approximation in the right hand side term is relevant at large distance (r >> tlc), 
where Vt*r » c. 

The relevance of the Alfven wings theory is based on a comparison of the wind's velocity, 
and the Alfven velocity. The Alfven speed in the wind's frame is 

V' A 2 = c- 2 + c A 2 = c- 2 + » p' /B> 2 , (4) 

where p' Q is the proper density of proper mass of the wind and B' Q is the unperturbed field 
observed in the wind's frame. The models considering a relativistic radial wind MHD 
model [Arons 2004, Kirk et al. 2009] imply an asymptotic Lorentz factor scaling as 

7oo ~ ^J /3 , (5) 
where o"o is the magnetization factor defined by 

a ° = — FT 6 

Then, it is shown in [MH] that the Alfvenic Mach number M' A = v r /V A scales as 1 — 
1/(2cJq^ 3 ), that is smaller than one. 

As the fast magnetosonic waves are even faster, the wind is slower than the fast magne- 
tosonic waves, and there is no MHD shock. Therefore, instead of being like a solar-system 
planet in the super- Alfvenic solar wind, the planets around a pulsar are more like Io in the 
sub-Alfvenic plasma flow of Jupiter's magnetosphere. Then, instead of being preceded 
by a bow shock, the pulsar's planet has a direct contact with the wind, and following 
the idea developed by Neubauer [1980], there is a pair of electric currents carried by two 
stationary Alfvenic structures called Alfven wings. They are attached to the planet on 
one side, and going far into the space plasma on the other side. The configuration of these 
electric currents is recalled in Fig. [TJ The engine of the Alfven wings is the convection 
electric field Eo = — vo x Bo induced by the motion of the wind into the magnetic field, 
that appears in the reference frame of the planet. It is generated by a potential drop U 
along the body of radius Rp, 

U = 2R P E = , (7) 

r 

where E = VqBq is directed perpendicularly to the wind flow and to the magnetic field. 

With each Alfven wing, there is an electric circuit consisting of a incoming current carried 
by the Alfven wing a coming from space, a current along the planet, and an outbound 
current going into space. The resistivity Saw associated to the Alfven wing was derived 
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from the MHD theory (see [MH]). It has a very simple form when the Alfven velocity is 
close to c, as in a pulsar's wind: 

Sa ~ — ■ (8) 
In that case, the order of magnitude of the total electric current is 



I- AW 



4(£ - Ei)R P E A , (9) 



where Rp is the planetary radius, and Ej, is an electric field along the planet caused by its 
ionosphere or surface finite conductivity [Neubauer 1980, MH]. According to Neubauer, 
the power dissipated by Joule effect along the ionosphere is maximized for matching of 
internal and external loads, that is, when E{ = Eq/2. In that case, for a P = 1 second 
pulsar and an Earth-like planet at 0.2 UA, I^w ~ 10 11 A. This current is not negligible. 
It has actually the same order of magnitude as the current that powers the whole pulsar 
magnetosphere [Goldreich and Julian 1969]. (Nevertheless, for asteroids, or for planets 
around recycled and less magnetized pulsars, the AW current is smaller than the Goldreich 
and Julian current by orders of magnitude.) 

The computation of this electric current involves a first invariant V s , derived in [MH]. 
With s = ±1, 

T 7 - sB(l-v -v/c 2 ) _ sB . , 

V s = V = V =; . (10) 

A(l - sB ■ v/c?) \^(l-sB -v /ci) 

The second equality is an estimate of the invariant taken from the unperturbed plasma 
area, at large distances from the wing. The parameter A can be expressed as a function 
of the unperturbed wind parameters, 

X=[ W ' + c^ l0 + c-^Bl ] 1/2 . (11) 

As seen with Eq. ([3]), at large distance from the light cylinder, B^ » B r and Vq = Vq t . 
Therefore, at the planetary distances, Bq ■ vq/c 2 << 1 and the invariant takes a simpler 
value, 

V s = v - ^. (12) 

The current given in Eq. ([9]) is parallel to this vector. Therefore, the angle of V s with 
the radial direction determines the geometry of the Alfven wing. The computation of the 
angle 9 between the wing and the radial direction is not done in [MH], and we present 
it in the following lines. Considering that at the planet's distance, the magnetic field is 
predominantly azimuthal and Vq ~ c, and 9 is given by the ratio of the two perpendicular 
components of V s , 

9 = arctan(£> /cA7o). (13) 

The magnetic field in A can be expressed as a function of the invariants of the wind, ^ 
and / given by the Eqs. and fl2]). In the algebraic development, / and ^ disappear, 
replaced by functions of the magnetization parameter <jq and the light cylinder radius 
r LC . Then, 



9 = arctan 



7o 



-2 \ 1/2' 



l + 0o7o +Vo(rLc/rf 



(14) 
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Considering the asymptotic value of the Lorentz factor given in Eq. (jSJ), and 7 >> 1, for 
a planet at large distance for the star, but still in the Poynting flux dominated wind, we 
find the very simple expression 

Ooo ~ 7» ~ ^o" 1/3 - (15) 

(In this last estimate, the effect of the residual magnetic field Bn, that is parallel to the 
flow, is completely neglected.) Figure [2] shows the possible aspect of the two Alfven wings 
connected to the planet, in the pulsar wind. The ambient magnetic field, that is spiral 
shaped (archimede spiral, with » B po i oidal ) is represented by field lines (thin lines on 
the figure). 

3 Radio emissions from the Alfven wings 

We can expect that with currents as strong as in Eq. (j^J), the Alfven wings associated to 
a planet in a pulsar's wind are strong radio-sources. 

These might be synchrotron radiation from the electrons that carry the current. But 
the current might be as well destabilized, and excited by transient Alfven waves caused 
by the fluctuations of the pulsar's wind passing along the planet. The conjunction of 
the destabilized current (filamentation, blobs with mirror points, density fluctuations) 
and transient Alfven waves might cause, for instance, acceleration and the production of 
cyclotron maser unstable distribution functions. Such coherent radiation process might 
be more efficient than the sum of the synchrotron radiation from all the particles. 

We have not studied yet the associated emission mechanisms that involve highly rel- 
ativistic plasmas; this topics is beyond the scope of the present report. But we can 
provide an estimate of critical frequencies that would play a role in such emission pro- 
cesses. Let us consider the planet "c" orbiting PSR 1257+12. The magnetic field, that 
is mainly azimuthal, can be estimated from Eq <^). The data in the tables allow to 
compute \& ~ 8.5 x 10 12 Wb, and the magnetic field at the distance a c of planet "c" 
is B (a c ) ~ 4 x 10~ 4 T. Therefore, the non-relativistic electron cyclotron frequency is 
f ce ~ 10 7 Hz. The synchrotron emission spectrum of a free electron, with a perpendicular 
component of the velocity corresponding to a Lorentz factor 7 e would increase up to the 
critical frequency f cr i t = fcell [Lang, 1999]. It is expected that the wind's electrons and 
the positrons created in the star's magnetosphere have already lost their perpendicular 
velocity, precisely through the synchrotron radiation. But the current could be carried, at 
least in parts, by particles extracted from the planet and accelerated. They could carry a 
non null perpendicular velocity. The spectrum of radiation of these particle would peak at 
f crit . Even for 7 e ~ 100 (that would represent a quite efficient acceleration), the spectrum 
would reach f cr u = 10 11 Hz: it would remains in the radio-frequencies range (decreasing 
beyond millimetre wavelengths). 
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4 Characteristics of the radio emissions from a planet around a 
pulsar 



A complete survey of the possible characteristics of the radio emissions of a planet around 
a pulsar remains to be done. Nevertheless, here are a few general remarks. 

Whatever their emission angle in the wind's frame we can predict their cone of emission for 
an observer in the frame of the Earth. As the sources are embedded in a ultra-relativistic 
wind, the emission angles would be contained in a cone of solid angle I/70 sr (relativistic 
aberration). 

It is important to notice that the emission would not come from the planet itself, but 
from the Alfven wings extending all along the planetary orbit. 

If we take into account the inclination of the star's magnetic field over its rotation 
axis, the magnetic field seen in the planet's orbital plane varies with the star's rota- 
tion angle. Therefore, the inclination of the wind's magnetic field angle with the orbital 
plane oscillate. As the Alfven wind makes a constant angle with the magnetic field, the 
Alfven wind inclination over the orbital plane oscillates too. The associated wavelength 
is A = ~ 2nc/Q*. Because of this oscillation, an observer would not see always 

the AW's radio emissions, they would appear pulsed. The observed period would be 

^observed ^* i ^orfr; (-^) 

where Q or b is the planet's orbital period (Q or b « ^*)- 

Another cause of variation of the period would be the Doppler shift associated to the 
varying distance of the sources (that are in the Alfven Wing) from the planet and from 
us. For emissions confined in the radial direction, this Doppler shift would be a function 
of the planet-star-observer angle. 

Because of the finite latitudinal extension of the AW (see Eqs. ffT4l) and (TT5]) ) and be- 
cause the radio-waves would be emitted in the radial direction, the possibility of their 
observation would depend on the observer-star-planet angle. As the latitudinal extension 
diminishes when the wind's Lorentz factor 7 increases, and because 70 is expected to 
be high, it is possible that the planetary radio emission are observed only during a small 
fraction of each orbit. The short laps of time of observability would be repeated regularly 
with the planet's orbital period. 



5 Conclusion 

A planet embedded in the highly relativistic, highly magnetized and sub-Alfvenic wind of 
a pulsar generates two pairs of electric current going inward and outward the planet. Such 
structures have been already observed in the case of Io, where the current flows between Io 
and Jupiter. These ribbons of current are the electric counterparts of stationary Alfven 
waves emitted by the planet. They are called Alfven wings. In the case of a planet 
in a pulsar's wind, the current is connected on one side to the planet, and fades into 
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Table 1: Data infered about the two stars known to host planets. 



Name 


P(s) 
Star's rotation 


ft* (s- 1 ) 


5* (Gauss) 
Magnetic 


M* (M ) 
Star's mass 


-R* (km) 
Star's radius 


PSR 1257+12 


0.006 


1010. 


8.8 x 10 8 


1.4 


~ 10. 


PSR 1620-26 


0.011 


567 


3. x 10 9 


1.35 


~ 10. 



Table 2: Data about the planets orbiting around pulsars. 



Name 


M P (M©) 


R P (i? ffi ) 


P orb (day) 


a (AU) 


e 


PSR 1257+12 a 


0.02 


0.28 


25. 


0.19 





PSR 1257+12 b 


4.3 


1.68 


66. 


0.36 


0.0186 


PSR 1257+12 c 


3.9 


1.62 


98. 


0.46 


0.0252 


PSR 1620-26 a 


794 


9,5 


36367. 


23. 





space under the control of the Alfven wave resistivity, that is of the order of the vacuum 
resistivity. 

The Alfven wing current amplitude has been estimated in [MH] . It is comparable to that of 
a pulsar (the Goldreich- Julian current) for a planet around a standard " 1 second" pulsar. 
It would be lower by three orders of magnitude in the case of recycled pulsars. Such 
systems of currents might be strong radio-waves emitter, especially for planets around 
young pulsars. The two Alfven wings would behave as two long antennas connected on 
one side to the planet. 

Because of the relativistic aberration, in the frame of an observer on Earth, the Alfven 
wings would emit radio waves mostly radially, in the outward direction (from the star to 
space). The emission angles would not describe the shape of a beam, but a fragment of 
a corona, beamed in the latitudinal direction, and in the azimuthal direction. The extent 
of the emission angle (6 in the latitudinal direction) would depend strongly on the wind's 
Lorentz factor 7 . 

Because of the (most expected) inclination of the pulsar's magnetic field over its rotation 
axis, the orientation of the magnetic field seen by the planet in the wind would be mod- 
ulated, with a frequency close to the neutron star's rotation period. As the orientation 
of the Alfven wings is determined by those of the magnetic field, their direction, and the 
direction of emission of the radio waves would wobble. This would result, for an observer 
in a fixed direction, into a pulsed emission. 

This study is preliminary. The mechanisms of the radio emission have not been studied, 
and a frequency spectrum is not proposed. We should just mention that this region is 
much less magnetized that the pulsar's inner magnetosphere where the "normal" radio 
emission are supposed to come from. Nevertheless, the physics of the radio emissions could 
not be transposed directly from Solar-system studies, because of the highly relativistic 
character of the pulsar's wind in which the sources would be embedded. 
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Table 3: Radii and corresponding keplerian orbital periods for the light cylinder of 1 s and a 10 
ms pulsars. 



Name 


P(s) 


r LC (km) 


Porb (s) 


Standard PSR 
Fast PSR 


1. 
0.010 


47 000 
477 


147 
0.147 




Figure 1: Schematic view of the induced electric field that generates Alfven wings. 
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